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ZJ-bursts from the Virgo cluster 
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Resonant annihilation of ultra-high energy cosmic neutrinos (UHECV) on the cosmic neutrino 
background (Cz/B) into Z bosons — the Z-burst mechanism — and its associated absorption and emis- 
sion phenomenology provide a unique, albeit indirect, probe of the Cz^B in its present state. In this 
paper, we examine the implications of gravitational clustering of the Cz^B in nearby galaxy clus- 
ters for the Z-burst phenomenology. In particular, we study the emission features of the Z-decay 
products originating from the Virgo cluster, and the potential of future cosmic ray experiments to 
observe clustering-enhanced Z-burst rates. We find that a detector with an exposure equivalent 
to three years of observations at the Extreme Universe Space Observatory (EUSO) will very likely 
measure these enhanced rates together with the associated UHECV flux, provided that the latter 
saturates current observational limits and the neutrino masses are quasi-degenerate, m Vi > 0.1 eV. 
In the case of UHECz/ fluxes below the electromagnetic cascade limit, or a hierarchical neutrino mass 
spectrum, an experimental sensitivity exceeding that of EUSO by at least two orders of magnitude 
is required to detect the clustering enhancements with any certainty. 



I. INTRODUCTION 



The existence today of a 1.95 K cosmic neutrino back- 
ground (Cz^B) — an exact analogue of the 2.73 K cosmic 
microwave background (CMB) — is a fundamental predic- 
tion of the standard big bang theory. Permeating the 
universe at an average number density of n Vi — n Pi ~ 
56 cm -3 per flavour, these relics of the big bang trace 
their origin to the freeze-out of the weak interaction when 
the universe was a mere one second old (T ~ 1 MeV), 
predating even the CMB photons by thirteen orders of 
magnitude in time. Yet, for the same reason that they 
decoupled so early, the Cz^B neutrinos have so far es- 
caped direct detection in a controlled laboratory setting. 
To date, cosmological measurements such as the CMB 
anisotropies and the large-scale matter power spectrum, 
and, independently, the observed light elemental abun- 
dances, provide the best probe of the C^B's presence in 
the early stages of cosmological evolution (e.g., pj). It 
is therefore our view that alternative avenues, however 
indirect, that could potentially afford us a glimpse of the 
Cz^B as it is today should be thoroughly explored. 

One such avenue is the Z-dip/burst mechanism and 
its associated phenomenology 0, which proposes to ex- 
ploit the resonant annihilation of hypothetical ultra-high 
energy cosmic neutrino (UHECz/) beams on the Cz^B as 
a target. Supposing such beams exist, the annihilation 
process zajheo + ^OB — > Z — > hadrons proceeds at the 
resonance energy E]^ with a cross section enhanced by 
several orders of magnitude compared to non-resonant 
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is a function of the neutrino and the Z masses, m Ui and 
mz, alone. The annihilation can be detected as absorp- 
tion dips in the incident UHECz; flux at E ~ ££ cs ("Z- 
dips") B0 and/or as emission features in the Z-decay 
products (nucleons and photons) 0. Indeed, in the lat- 
ter case, the happy coincidence between E™ s ~ 10 21 eV 
(for m Vi ~ 1 eV) and the energies of the most energetic 
cosmic rays observed by AGASA @, Fly's Eye @, HiRes 
Q, Yakutsk @, and Haverah Park @ has long led to a 
possible identification of ultra-high energy cosmic rays 
(UHECR) above the Greisen-Zatsepin-Kuzmin (GZK) 
cut-off energy Eqzk ~ 4 x 10 19 eV with Z-burst 
nucleons and photons 0, ^| . 

Clearly, the success of the Z-burst mechanism as a 
means to detect the Czd3 and/or as an explanation of 
UHECR depends first and foremost on the UHECz; fluxes 
(the beam); Figures ^ and [21 summarise the current sta- 
tus of the search for such fluxes, the projected sensitiv- 
ities of ongoing and planned experiments, and predic- 
tions from various theoretical models. A second factor 
is the nature of the Czd3 density distribution (the tar- 
get). To this end, it is important to note that an oscil- 
lation interpretation of the atmospheric and solar neu- 
trino data (e.g., [2^) implies that at least two of the 
mass eigenstates in the Cz^B are nonrelativistic today, 
i.e., m Vi > T v ~ (4/H) 1 / 3 ^ - 2 K - 10~ 4 eV, and at 
least one mass must exceed \J Am 2 tm ~ 0.05 eV. These 
nonrelativistic neutrinos are subject to gravitational clus- 
tering on existing cold dark matter (CDM) and baryonic 
structures, possibly causing the local CzTJ density to de- 
part from the cosmological average. 

Standard large-scale structure theories tell us that, in 
the currently favoured ACDM cosmology, {O m , SI a, h} = 
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FIG. 1: Upper limits on the diffuse neutrino flux per flavour 
F Va + F Da ,a = e,n,r, from RICE GLUE [3, Fly's 
Eye+AGASA Q, FORTE IT|] and AMANDA (assuming an 
E~ 2 differential spectrum) [lql . assuming full flavour mixing 
en route to Earth. Also shown are the projected sensitivities 
of Auger in v e , v u modes and in v T mode (bottom swath) 
E3, ANITA Jig, EUSO E3, IceCube and SalSA [H, 
i.e., one event per energy decade for the indicated duration. 
A prototype of the full ANITA experiment, ANITA-lite [2^1 . 
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{0.3,0.7,0.7}, fluctuations in the CVB ought to track 
their CDM counterparts at scales above the neutrino 
free-streaming length (e.g., 0). The inferred local large- 
scale matter distribution from peculiar velocity measure- 
ments p9| therefore precludes any increase in the relic 
neutrino content of the local GZK zone (~ 50 Mpc) by 
more than a factor of two due to gravitational cluster- 
ing alone [3(j. However, clustering in local gravitational 
potential wells such as galaxies and galaxy clusters may 
still be sizeable on the sub-Mpc to Mpc scale. Indeed, for 
neutrino masses satisfying experimental and cosmological 
bounds, 1 detailed calculations show that the CvB over- 
densities in the largest galaxy clusters (~ 10 15 M Q ) can 



1 Limits on the neutrino mass derived from cosmological measure- 
ments, especially the matter power spectrum P(k), depend on 
the data sets used and the priors assumed. For example, analy- 
ses including also the Lyman a data and information on the dark 
matter-galaxy bias tend to produce bounds that are much tighter 
than those derived from the shape of P(k) alone. Furthermore, 
the neutrino mass appears in combination with several other un- 
certain cosmological parameters in the determination of P(k), 
such as a running scalar spectral index n s (k) and the effective 
number of thermalised fermionic degrees of freedom N v . These 
degeneracies can lead to considerable relaxation in the bound on 
m Ui . In the present wor k, w e take a conservative upper limit of 
^jiHi/j < 1.8 eV (2cr) I'S ll for three degen erat e neutrinos, de- 
rived from the SDSS galaxy power spectrum 1321 and the WMAP 
data l33l assuming a constant n 3 and N v = 3. La bor atory 
bounds from tritium /3-decay experiments, m„ < 2.2 eV l34l , and 
from neutrinoless double beta decay, m„ < (066-^2.70) eV p|. 
are not yet competitive. See the reviews [il l 3(1 



be as much as a thousand within the central ~ 100 kpc re- 
gion [23, H3 • An immediate consequence for the Z-burst 
scenario is a possible directional dependence in its emis- 
sion features, even if the UHECV sources are isotropically 
distributed; the highest number of Z-burst events should 
originate from the directions of nearby galaxy clusters. 

The implications of gravitational neutrino clustering 
for the Z-burst emission spectra were first investigated 
quantitatively in reference |37| . In the present paper, we 
extend the said analysis in several ways: 

1. We present a more accurate determination of the 
Cz^B overdensities in and around galaxies and 
galaxy clusters based on the calculations of refer- 
ence 30] . which take into account nonlinear effects 
in the clustering process. The linearised method 
adopted in (3?J systematically underestimates the 
overdensities by a factor of several particularly in 
the virialised region of large galaxy clusters. 

2. We compute both the primary nucleon and pho- 
ton spectra. A generic feature of the Z-burst sce- 
nario is a photon to nucleon ratio of almost 20:1 
at the Z production site. This ratio is greatly re- 
duced upon arrival at Earth because of the much 
shorter attenuation length for photons than for nu- 
cleons. Still, unless the photons are heavily at- 
tenuated by a strong universal radio background 
(URB) and/or by strong extragalactic magnetic 
fields ( > 10~ 10 G), we expect a predominance of 
photons in the observed primaries. 



3. As in [37|. we assess the experimental prospects 
for observing enhanced Z-burst rates from nearby 
galaxy clusters, focussing on our nearest neighbour: 
Virgo. At an average distance of ~ 15 Mpc from the 
Milky Way and with a mass close to 10 15 M Q pj. 
the Virgo cluster is able to accumulate a sizeable 
excess of relic neutrinos while sitting close enough 
such that the Z-burst primaries suffer minimal en- 
ergy loss during their propagation to Earth. If 
clustering-enhanced Z-bursts were to be seen at all, 
the Virgo cluster would be the prime site. 

For rate estimates, we shall take the Extreme Universe 
Space Observatory (EUSO) experiment [23] as our fidu- 
cial detector. The sensitivity of EUSO to UHECR and 
UHEC^ above 10 20 eV is about three orders of magni- 
tude beyond what is available to date from AGASA and 
HiRes. As presently designed, EUSO will consist of a 
two-metre Fresnel lens positioned on the International 
Space Station at a height of 400 km. The lens will focus 
near-UV fluorescence emitted by radial de-excitation of 
N2 in the air shower. EUSO operates as a space-based 
eye looking for tracks in one gigantic cloud chamber: the 
Earth's atmosphere. The lens' opening angle will be 60°, 
giving the instrument an enormous field of view (FOV), 
7T x (tan 30° x 400 km) 2 = 1.7 x 10 5 km 2 . Since the at- 
mospheric density decreases exponentially with altitude 
with an 8 km scale-height, the EUSO FOV encompasses 
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the task at hand. 

The paper is structured as follows. In section[n]we de- 
scribe the calculational procedure for the Z-burst fluxes. 
Section II I II discusses the modelling of the C^B density 
distribution. Our predictions for the Z-burst fluxes orig- 
inating from the Virgo cluster and the corresponding 
event rates expected at EUSO are presented in sections 
FVl and W\ respectively. We conclude in section IVTl 



II. PRELIMINARIES 



The differential flux in particle species ip (ip = j,p+n), 



FIG. 2: Prospects for ultra-high energy neutrino detection 
in the next decade (fluxes required to produce one event per 
energy decade). For the year 2008 (long dashed/red line), we 
assume three years of Auger data and one 15-day ANITA 
flight. For 2013 (dash-dot/blue line), we assume 8/3/3/4 
yr Auger/EUSO/IceCube/SalSA, and three ANITA flights. 
These curves are indicative only. The sensitivity will improve 
if we consider also the Westerbork radio observatory |23f . if 
further projects such as Auger North and OWL [24| are re- 
alised, or if the EUSO or ANITA flight times are extended. 
Also shown are a wide sample of UHEC^ flux predictions, the 
current observational upper bound (solid shade) from Fig- 
ure Q and the electromagnetic cascade limit from EGRET 
[SI |2(|- The points labelled "Z-burst" (green) denote the 
UHECf fluxes required to explain the post-GZK events ob- 
served by AG AS A in terms of Z-burst secondaries |2^| . 



an equivalent of 1.4 x 10 6 km 3 of air at surface density, 
corresponding to more than a Teraton of mass. 

Due to the thin height of the Earth's atmosphere, the 
distance d between the EUSO lens and the observable 
air shower is highly constrained. For a shower directly 
below the lens, d — 400 km to within a few percent; for 
an event on the periphery of the FOV, d increases by only 
l/cos30° = 1.15. Thus, the energy threshold for EUSO, 
determined by the 1/d 2 fall-off of signal, is quite sharp. 
In its presently planned configuration, the threshold of 
EUSO is a few times 10 19 eV, fortuitously positioned to 
observe the region at and above the GZK cut-off. 

The duration of the EUSO experiment is nominally 
three years, but extensions beyond that seem likely. On 
the other hand, the smaller Auger Project |40J will likely 
collect data for a decade, so that its exposure (product of 
acceptance and time) may come to rival that of EUSO. 
Unfortunately, the present configuration for Auger is lim- 
ited to the Southern Hemisphere, casting a blind eye in 
the direction of Virgo. Other existing cosmic ray facilities 
are not competitive with the size of the EUSO FOV. The 
angular resolution expected for EUSO showers is about 
one degree, or 20 mrad. Other neutrino telescopes also 
offer similar resolutions. As we shall show later, the CVB 
"halo" expected for the Virgo cluster spans several to ten 
degrees, so a one-degree resolution is sufficiently fine for 



Fil>\z{E, 9, 



dN, 



tp\Z 



dE dA dt dn ' 



(2) 



denotes the number of ip particles arriving at Earth with 
energy E in the direction {9, </>}, where {0,0} labels the 
centre of the Virgo cluster (i.e., galaxy M87), per unit 
energy per unit area A per unit time t per unit solid 
angle Q. In the Z-burst scenario, this is given by [27j 
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where F Vi (E Vi , r) is the ith mass UHEC^ flux 2 at energy 
E Vi at the Z production point at a "look-back" time 
t = r/c, n Vi (r,9,4>) the Cz/B number density, c^p^s) 
the Z production cross section at centre-of-mass energy 
sfs = ^/2m v% E v%1 Br(Z -> hadrons) = (69.89 ± 0.07) % 
the branching ratio, dN^/dE^ the energy distribution of 
the produced tjj particles with energy E^, , and the propa- 
gation function (r, E^ ; E) gives the expected number 
of ip arriving at Earth with energies above the threshold 
E per particle created at r with energy E^,. 3 



2 Technically, one should write F Vi as p \ Upi \ 2 , where 
|t^3;| 2 = I < va\vi > 1 2 is the projection probability for flavour 
eigenstate /3 onto mass eigenstate i, and F V/j is the initial flux 
of flavour [3. However, this exactness introduces an unnecessary 
layer of detail for the investigation performed here. 

3 In the process of becoming nonrelativistic, the CVB neutrinos 
lose their handedness and depolarise to populate all of their spin 
states. If neutrinos are Majorana particles, as is predicted by the 
most neutrino mass models, then all CVB spin states will par- 
ticipate in annihilation, and equation is correct as written. 
However, if neutrinos are Dirac particles, then half of the nonrel- 
ativistic spin states are "sterile" states, and the rate presented 
in should be reduced by a half. 
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Since the cross section a I/i p i (s) is sharply peaked at the 
resonance energy s — m 2 z , we may approximate 

dE Vi F Ui (E Ui )a ViPi (a) ~ E™F Ui (Elf) (a ann ) , (4) 



where the superscript "res" denotes resonance, and 
(cann) = J ds a ann /m z — 40.4 nb is the energy-averaged 
s-channel Z-exchange annihilation cross section. The 
UHECV fluxes are modelled as 



F Vi (E Ui ,r) = F Vi {E Vii 0)[l + z(r)] a , 



(5) 



given our insufficient knowledge about their sources. 
Here, F Vi (E Vi , 0) is the neutrino flux incident on Earth, 
the index a characterises the sources' cosmological evo- 
lution, and the redshift z(r) is related to the look- 
back time via dz = —(1 + z)H(z) dr/c, with H(z) — 
Hq [f2 m (l + z) 3 + SIa] 1 / 2 linking the Hubble parameter 
at z to its present value, Hq = h 100 km s _1 Mpc -1 . Fur- 
thermore, we take n Vi (r, 9, <p) ~ n Pi (r, 9, </>), 4 since signif- 
icant relic neutrino-antineutrino asymmetries ( > 1) are 
incompatible with big bang nucleosynthesis in the pres- 
ence of bi-large mixing inferred from the solar and atmo- 
spheric neutrino experiments |4lj|. 5 Thus, together with 
the definition F™ s = F Vi (E™ s , 0) + F Pi (E™ s , 0), equation 
© now becomes 

F^ z (E,6,ct>) = ^ 2 Br(Z ^ hadrons) (<r anll ) F™ s 

i 

x dE dr (l + z) a n Vi {r,6,<t>) 
Jo Jo 

dP^E^-E) 



dE 



(6) 



The functions Q-d>(y) — E v /2 • dN^,/dE^, with y — 
Am l/ E^/m 2 z , are the boosted momentum distributions 
from hadronic Z-decay, normalised to (N p+n ) = 2.04 for 
ijj = p + n, and to (N 7 ) = 2(N„o) + (N„±) = 37 for 
ip = 7. 6 Detailed forms for (y) can be found in refer- 
ence |27j . 



4 It is interesting that the presence of a neutrino-antineutrino 
asymmetry increases both n l>i + na i and \n v , — noA, while at 
the same time driving one of n Vi or n„ i to zero exponentially. 
This has the curious consequence of increasing the Z-burst rate, 
but not necessarily the Z-dip depth. For the latter, the severe 
suppression of n Vi or na i in the CVB gives an asymptotic depth of 
one half for the Z-dip. Of course, the increased number of events 
concomitant with a neutrino-antineutrino asymmetry improves 
the statistics of the dip. 

5 Constraints on the cosmological and v T neutrino-antineutrino 
asymmetries are applicable insofar as large-angle u e <-> i>^ lT oscil- 
lations are operational prior to neutrino decoupling. An obvious 
way to evade these bounds is to suppress these oscillations by 
way of new, non -standard matter effects. One such scenario is 
presented in |42| . in which the suppression arises from a hypo- 
thetical flavour-dependent neutrino-majoron coupling. 

6 The photon count includes also electrons and positrons from 
charged pion decay. These are relevant for the development of 
electromagnetic cascades. 



The propagation functions P^ (r, ; E) account for 
the interactions encountered by ip between its production 
point and Earth. For ip = p+n, energy loss arises primar- 
ily from pion and e + e~ production through nucleon scat- 
tering on the CMB. The corresponding P p+n (r, E p+n \ E) 
has been calculated in detail in ji^ \ an d is publicly avail- 
able at For the computation of P~(r,E~;E), we 
adopt the continuous energy loss approximation, which 
asserts that the photon energy degradation proceeds as 



dE 
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dr 
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e z (E) i + z 
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in which l z (E) = (1 + z)- 3 £ (E(l + z)), and £ {E) is 
the photon energy attenuation length due to pair and 
double-pair production on the diffuse extragalactic pho- 
ton background, and inverse Compton scattering of the 
produced pairs. Values for £q(E) incorporating various 
ass ump tions about the poorly known URB can be found 
in |45j and are summarised in Furthermore, we 

assume the number of photons N 7 to be constant at 
energies > 10 18 eV due to the small inelasticities in 
this energy range. Below ~ 10 18 eV, N 7 increases as 
dN-y — N 7 dr /£ Z (E), so as to maintain energy conserva- 
tion (excluding losses due to the universal expansion). 



III. NEUTRINO DENSITY DISTRIBUTION 

The Cz^B density distribution in the universe is taken 
to be uniform and at the cosmological average n Vi = 
n ei (~ 56 cm -3 at z = Q), except in the vicinity of a 
galaxy cluster, in which case we use the neutrino number 
densities provided in reference |30| (reproduced here in 
Figure These densities are obtained by solving, with 
a particle realisation, the Vlasov and Poisson equations, 



dr 



P 



df 

am Vi dx 

72 J. A _/"r„2 



am Vi \I<f> 



0, 



AttGci [p m (x,T) 



dp 

- Pm{T)] 



(8) 
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Here, f(x,p,r) is the neutrino phase space distribu- 
tion, where x, p, and r are the usual comoving coor- 
dinates, conjugate momentum, and conformal time, re- 
spectively; a is the scale factor, p m the mean universal 
matter density, and we assume a ACDM cosmology with 
{fl m ,n A ,h} = {0.3,0.7,0.7}. The halo density profile is 
taken to be of the Navarro-Frenk- White form [46J , 
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where r is the radial distance from the halo centre, and 
the parameters r s and p s are determined by the halo's 
virial mass M v ; r and concentration c via 
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FIG. 3: Relic neutrino number density per mass eigenstate, function of the radial distance from the halo centre, 

for halo virial masses (top to bottom) 10 15 , 10 14 , 10 13 , 10 12 , in units of (0.7 /h) M@, and neutrino masses indicated in the 
figure. All curves are normalised to the expected mean density n Ui = n„ i ~ 56 cm -3 . 



In addition, high resolution simulations of reference [47 
provide a tight correlation between c and M V i T , 



c(z) 



1 + z \1.5x 10 13 h- 1 M Q 



-0.13 



(13) 



The net result is a dependence of the halo density profile 
on the halo's virial mass alone. 



IV. Z-BURST FLUXES FROM VIRGO 

To estimate the Z-burst fluxes from the Virgo clus- 
ter, we assume the cluster's virial mass to be 10 15 Mq, 
centred on galaxy M87 at a distance D ~ 15 Mpc from 
Earth. We integrate equation © along the line of sight 
in the direction {9, </>}. The upper integration limit is 
taken to be i? ma x = 3000 Mpc, although the choice of 
Rmax has little impact on the results, provided R max ex- 
ceeds the GZK distance ~ 50 Mpc. We note that contri- 
butions from gravitational clustering in the Milky Way 
halo (~ 10 12 AT Q ) are generally negligible despite its prox- 
imity. This is because enhancements in the C^B density 
therein are no more than a factor of twenty even for the 
most massive neutrino considered here, and regions of 
substantial overdensity (n Ui /n Ui > 2) are limited in ex- 
tent (< 100 kpc). 

Since the CvQ distribution is not uniform inside the 
Virgo cluster but decreases with the halo radius, we ex- 
pect the Z-burst fluxes to vary with angular distance 9 
(recall that {9 = 0, <f> = 0} labels the centre of the clus- 
ter) in the same manner. Figure 0] shows the predicted 
nucleon and photon fluxes as functions of E and 9 for a 
range of neutrino masses and three different URB scenar- 
ios documented in the literature. The choice of a (i.e., 



the source evolution parameter) is irrelevant for Virgo 
Z-bursts, since the Virgo-Earth distance corresponds to 
a mere z ~ 0.003: we take a = 0. 



For definiteness, we have opted to evaluate the UHECf 
fluxes F™ s at the electromagnetic (EM) cascade limit, 
although fluxes larger by a factor of five to twenty-five, 
depending on the neutrino mass, are permitted by obser- 
vations. The EM cascade limit applies generally to trans- 
parent sources wherein neutrinos are produced in a chain 
of particle decays (e.g., pions, W, and Z). The decay 
process inevitably generates also photons and/or elec- 
trons with an energy fluence comparable to that of the 
neutrinos. Subsequent EM cascades via collisions with 
the diffuse extragalactic photon background, notably the 
CMB, bring the cascade photons into the energy range 
30 MeV-=- 100 GeV probed by EGRET [H . Observation 
of the diffuse extragalactic 7-ray background at these en- 
ergies therefore places a convenient upper bound on the 
diffuse UHEC^ fluxes H|. This is known as the EM 
cascade limit (or cascade limit, for short). 7 



For comparison purposes, we include in Figure 0] also 
the spectrum of cosmic ray protons produced at diffuse 



7 The EM cascade limit presented in [26 1 (displayed here in Figure 
Ul derives from the estimate of the diffuse 7-ray background pub- 
lished in l25l . Recent works |411 suggest that the latter analysis 
may have overestimated the extragalactic contribution to the 7- 
ray background by roughly a factor of two. The EM cascade 
limit may therefore be stronger correspondingly. 
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FIG. 4: Z-burst nucleon and photon fluxes from the Virgo cluster {M v i t = lO 15 A/0, D = 15 Mpc) per neutrino species for three 
URB scenarios (minimal, moderate, and high), assuming the EM cascade limited UHECf fluxes. Red (solid) lines represent 
photon fluxes originating from points located at angular distances (top to bottom) 9 = 0°, 4°, 10°, and 180° (i.e., an unclustered 
C^B) from galaxy M87. Blue (dashed) lines denote the corresponding nucleon fluxes. The Akeno+AGASA (upper) and Fly's 
Eye+HiRes (lower) normalised extragalactic proton fluxes are shown in green (dotted lines). 



extragalactic sources, 



F plbkd (E;A,f3) = / dE, 



dP p (r,E p ;E) 



dE 



(14) 



where {A,n, 7} = {2.37 x KT 10 eV^m^s- 1 , 3.65, 2.54} 
and {1.25 x KT 10 eV^m^s- 1 , 3.45, 2.54} are the best 
fit parameters for the existing Akeno+AGASA and 



Fly's Eye+HiRes data, respectively, in the energy range 
10 17 ' 6 H-10 20 eV :50|. A lower integration limit of R min = 
50 Mpc has been imposed in the evaluation of F p \ hkd , 
since no known sources reside within this distance. The 
resulting spectrum exhibits an accumulation at the GZK 
scale 4 x 10 19 eV, and a sharp drop beyond. 

The GZK suppression of nucleons from distant sources 
provides a clean environment for the study of Z-burst 
events originated by cosmic neutrino messengers. For 
the cascade limited UHECV fluxes assumed here, it is 
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clear from Figure 0] that Z-burst events dominate above 
10 20 eV, especially in the direction of Virgo. Note that 
substantial enhancements in the Z-burst fluxes due to 
CVB clustering on Virgo can still be seen at ~ f 0°, well 
beyond the cluster's visible region < 5°. This follows 
from the assumption of an extended CDM halo, which 
accumulates neutrinos in the outer region gravitationally 
despite its seeming invisibility. 

In Figure we present the "sky map" of the energy- 
integrated Z-burst fluxes for nucleons plus photons above 
E t h = 2 x 10 20 eV, centred on the Virgo cluster: 

dEF^ z (E,6). (15) 

The shape of this curve is not very sensitive to the choice 
of the URB, but is highly dependent on the neutrino 
mass, since the latter determines the amount of gravita- 
tional clustering available to the CvQ. In principle, the 
angular distribution and the energy dependence of the 
Z-burst events provide independent confirmation of the 
neutrino mass (and a consistency check of the Z-burst 
mechanism). The angular resolution of EUSO and other 
UHECR detectors is typically one degree, sufficient to 
map out the shape of the extended Virgo halo. 



V. EXPECTED EVENTS AT EUSO 

The EUSO experiment detects UHECR by observing 
near-UV fluorescence emitted by nitrogen molecules in 
the extensive air showers generated by the primary cos- 
mic ray particles' interactions with the Earth's atmo- 
sphere. For primary nucleons and photons, the cross 
section for interaction is sufficiently large that the inter- 
action takes place high in the atmosphere with near unit 
probability. Thus, the acceptance of EUSO for nucleons 
and photons, A, is just the projected FOV normal to the 
direction of the source, i.e., \x FOV ~ 0.85 x 10 5 km 2 , 
times the solid angle on the sky of the emitting re- 
gion, times another factor ^ accounting for the block- 
age of the upcoming beam by the opaque Earth, times 
the duty cycle of the instrument, i.e., the fraction of 
time "on" . The EUSO instrument can record fluores- 
cence signals only on moonless nights devoid of high cir- 
rus clouds. The duty cycle for such clarity is estimated 
to be 20 %. This amounts to an all-sky acceptance of 
-44tt = 1.1 X 10 5 km 2 sr, including the duty cycle fac- 
tor. The nominal duration of the EUSO experiment T 
is three years (~ 10 8 s), so that an all-sky exposure of 
Sin = A^T ~ 1.1 x 10 13 km 2 s sr is anticipated. If 
the lifetime of EUSO is extended, the exposure will be 
proportionately larger. 

For the Virgo cluster, the expected number of Z-burst 
primaries in the energy interval (Ej, -Ey+i) originating 
within an angular distance 9 from the cluster centre is 




0.8 r 
0.6 - 
0.4 r 



0.2 r m,=0.07 eV j: 

t_i i i I i i i I i i i I i i i I i i i ~. 

2 4 6 8 10 

6 [degrees] 

FIG. 5: Energy-integrated Z-burst fluxes (nucleons plus pho- 
tons) above — 2 x 10 20 eV as functions of angular distance 
from the centre of Virgo, normalised to their corresponding 
values at 9 — 0°. Red (solid), blue (dashed) and green (dot- 
ted) lines denote, respectively, the minimal, the moderate and 
the high URB scenarios. 



given by 



%W = T^/ 3+1(1E l d^' ! dO'sme'F^ z {E,0',<l>'). 
4tt J e . J J 

(16) 

Figures El and [7] show the integral l|16|l evaluated for 
6 = 10°, assuming, respectively, UHECz/ fluxes at the 
current observational and the EM cascade limits, with 
j = 1,...,4 designating the four logarithmic bins be- 
tween 10 20 -^10 22 eV. Also displayed in the figures are the 
numbers of extragalactic background protons anticipated 
in the same solid angle Ail = 27r(l — cos 9) ~ 0.1 ster, 
and of Z-burst events in the same Ail for an unclustered 
Cj/B. 
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FIG. 6: Number of Z-burst events per neutrino species originating within 9 — 10° from M87 expected at EUSO in three 
years, assuming UHECV fluxes at the current observational limits. Solid red squares denote the photon events, while solid blue 
triangles refer to the nucleons. The corresponding predictions for an unclustered C^B in the same solid angle ASl ~ 0.1 ster 
are indicated by the open red squares and open blue triangles. The number of extragalactic background protons anticipated 
in Afl, assuming the Akeno+AGASA (upper) and the Fly's Eye+HiRes (lower) normalisations, are represented by the green 
crosses. 



A. Most optimistic scenario: hidden sources 



quasi-degenerate, m Vi > 0.1 eV. 8 The resulting number 
of events at EUSO in the ~ 0.1 ster solid angle can be 



As is evident in Figure El the most favourable circum- 
stances under which clustering-enhanced Z-burst emis- 
sions from the Virgo cluster may be observed occur when 
(i) the UHECV fluxes saturate the current observational 
bounds, and (ii) the neutrino masses are "large" and 



The mass splittings inferred from the atmospheric and solar neu- 
trino oscillation experiments |28| . Am^ tm ~ 2 X 10 — 3 eV and 
Am 2 un ~ 7 X 10 — 5 eV 2 , respectively, imply that the three mass 

eigenstates are quasi-degenerate when m„ S> u Am 2 tm . 
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quite large, totalling ~ 100 events in the energy decade 
10 20 10 21 eV for m Vi = 0.6 eV, ~ 40 for m Vi = 0.3 eV, 
and ~ 10 for m Vj = 0.15 eV, for three degenerate species. 
The corresponding numbers for an unclustered CVB are 
~ 6, ~ 4, and - 2. 

Theoretical models capable of generating neutrino 
fluxes at these formidable energies and magnitudes — 
either by way of "bottom- up" astrophysical accelerators, 
or through "top-down" decays of super heavy particles 
from beyond the standard model — are not lacking (e.g., 
|5l)'). The significant hurdle facing these UHECz^ sources 
is that they must be opaque to nucleons and high energy 
( > 100 MeV) photons, in order not to exceed the diffuse 
7-ray background observed by the EGRET experiment 
(cf. the EM cascade limit described in section Hv|l . A 
proof-of-principle candidate is a mirror (hidden) topolog- 
ical defect coupled to the standard model (SM) through 
mirror-SM neutrino oscillations; by construction, the re- 
sulting SM neutrino fluxes is free to saturate the upper 
observational limits Furthermore, because photons 
are regenerated by Z-decays, their subsequent EM cas- 
cade down to the energy range probed by EGRET can 
lead to additional constraints on the UHECj/ sources, 
particularly on the sources' cosmological evolution pa- 
rameter a. For UHECi' fluxes saturating current ob- 
servational limits, sources with a < are consistent with 
the EGRET bound 53] . 9 The next generation of UHECz/ 
observatories will put these neutrino source scenarios to 
a definitive test (cf. Figure^- 

If neutrino masses are not quasi-degenerate, but rather 
hierarchical with m Vi < 0.1 eV, then their clustering in 
Virgo is small or negligible. For m Vi = 0.07 eV, the 
clustered and the unclustered rates differ generally by no 
more than a factor of two (cf. Figures [S] and HjJ . Judging 
from our numbers in Figure (bottom row) , even for 
an optimised UHECf flux and the best URB scenario, 
an experiment with at least ten times the exposure of 
EUSO (three years) is required to record one clustering- 
enhanced Z-burst event from Virgo, another two times 
to record an event within Af2 with no CiaB clustering, 
and another five times to resolve the difference between 
the two at lcr. 

Of course, if there is indeed a large UHECi/ flux, it can 
also be measured directly by EUSO. It is a simple matter 
to obtain a good estimate of the EUSO acceptance and 
exposure for neutrinos from those for cosmic ray nucleons 
and photons. To a very good approximation, the atmo- 



The EGRET constraint on a in '531 applies only to a homoge- 
neous CVB. However, we do not expect it to be too seriously 
affected by enhanced Z-burst emissions from gravitational clus- 
tering at Virgo. This is because the enhancements are no more 
than a factor of twenty, and are concentrated in a small solid 
angle, Af2 ~ 0.1 ster, less than one hundredth of the whole sky. 
Changes in the all-sky emission rate due to enhanced emissions 
from Virgo can only be at the < 15 % level. Other galaxy clus- 
ters of comparable sizes are likely too far (i.e., > 50 Mpc, outside 
the GZK zone) to be of great concern. 



spheric density decreases exponentially vertically, with a 
scale-height of xt ~ 8 km. Thus, the effective volume is 
XhX FOV. The nucleon density in this volume is p{x = 
0) x Nai where p(x = 0) is the atmospheric density at 
sea level in g cm -3 , and Na is Avogadro's number. Thus, 
relative to the aperture for cosmic rays, we have for neu- 
trinos an extra factor of 2xh p{0) Na <j ,,n, where <j u n rs 
0.77 x 10~ 31 (E u /10 20 eV) 36 cm 2 Q is the neutrino- 
nucleon cross section, and the factor of two arises be- 
cause for each neutrino with an oblique trajectory, the 
extra path length compensates the reduced projection of 
the FOV. Putting in numbers (e.g., the vertical slant- 
depth, Xh p(0) — 1030 g cm~ 2 , is a well known number), 
one arrives at A% = 0.9 x 10~ 4 (E„/10 20 eV) - 36 » 
10~ 4 x Aatt at E v = 10 20 eV, with an increase in A\ v of 
2.3 per decade of energy beyond 10 20 eV. The estimate of 
this factor assumes that the experimental efficiencies for 
cosmic rays and for neutrinos are the same. 10 For EUSO, 
the efficiencies for cosmic rays and for neutrinos are each 
near unity at energies above 10 20 eV. The discriminator 
between cosmic ray and neutrino initiated events is the 
depth of the origin of the shower in the atmosphere. For 
hadrons and photons, the shower begins high in the atmo- 
sphere, while for neutrinos it begins much lower, where 
the air is densest. 

In its all-sky search, EUSO should record some ~ 
E v F V {E V ) UHEC^ events at energies above E v . In 
the energy interval 10 21 10 22 eV, this is of order 200 
events per neutrino flavour for F v at the observational 
limit, and ten times fewer for F v at the cascade limit. 
This measurement will establish the absolute normali- 
sation of the UHEC;/ flux, thereby removing the final 
uncertainty in the Z-burst calculation. It is therefore 
probable that clustering-enhanced Z-burst emissions and 
UHECi^ events will be simultaneously measured as soon 
as EUSO has completed its nominal three- year flight, if 
conditions (i) and (ii) are indeed satisfied. In such a case, 
the associated Z-dips [5f| in the incident UHEC^ flux will 
also have been resolved by the next generation of dedi- 
cated UHEC^ detectors such as ANITA in the same time 
frame. Remarkably, viable UHECz^ fluxes and neutrino 
masses in the ~ 0.1 eV ballpark are also able to produce 
the correct amount of Z-burst nucleons and photons to 
explain the cosmic ray events observed by AGASA at 
energies above the GZK cut-off. 

It is interesting to note that if EUSO, or any succes- 
sor experiment, succeeds in measuring Z-bursts and the 
diffuse UHEC^ rate at E r ° s : then the neutrino-nucleon 
cross section o v n can be inferred at E r J: s ~ 10 22 eV, far 
above energies available to terrestrial accelerators. The 
ratio of the Z-burst rate and the diffuse UHEC^ rate 



This same factor describes the relative acceptances and exposures 
of any experiment triggering on downcoming and horizontal at- 
mospheric events. Besides EUSO, another large area example is 
the Pierre Auger Observatory. 
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FIG. 7: Same as Figure El but for the EM cascade limited UHECV fluxes. 



is independent of the flux, but dependent on the two 
cross sections, (<7 ann ) and <j v n- The former is determined 
purely by weak interaction physics to be 40.4 nb, leaving 
only ov at, which depends on the QCD dynamics of the 
nucleon, as the unknown variable. 



B. Less optimistic scenario: transparent sources 

For transparent sources, the EM cascade limit on the 
UHECV fluxes applies, thereby excluding Z-burst emis- 
sions as an explanation of the AGASA post-GZK excess 
(cf. Figure [2J|. Nevertheless, EUSO has a projected ex- 



posure three orders of magnitude larger than the expo- 
sures of existing experiments; one may therefore hope 
for the eventual discovery of Z-burst events in the fu- 
ture, even if Z-bursts are not the source of the AGASA 
events. However, with a cascade limited UHECj/ flux, 
the observation of Z-burst rates in the direction of Virgo 
seems difficult to realise even with three years of EUSO. 
This remains true even with the substantial clustering en- 
hancement available to quasi-degenerate neutrinos (Fig- 
ure UJ. Better statistics may be achievable by widening 
the solid angle about Virgo, but inevitably at the ex- 
pense of narrowing the tell-tale gap between the clustered 
and unclustered rates. Furthermore, to pursue this low 
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signal-to-noise scenario, one would likely require a more 
accurate modelling of gravitational neutrino clustering in 
the local universe — something beyond our simple halo de- 
scription. All in all, a larger experiment is required. The 
OWL and multi-OWL proposals _24j would put (multi) 
satellites into orbit, to provide a FOV that dwarfs even 
EUSO. Experiments using LOw- Frequency radio antenna 
ARrays (LOFARs) have also been proposed. The proto- 
type for this kind of detector has just reported positive 
identification of UHECR H|. The Westerbork [23| LO- 
FAR facility in the Netherlands may offer the potential 
to improve the event rates of UHECR and UHEC^ by 
two or more orders of magnitude beyond EUSO. 

Suppose the UHECi^ flux comes from transparent 
sources, in which case it is cascade limited and cannot 
explain the post-GZK events in AGASA. Further sup- 
pose that the post-GZK AGASA spectrum and rate are 
correct. Then, is the "background- free" window for neu- 
trino physics above Eqzk closed? The answer is, not 
necessarily, if the Cz/B does cluster sufficiently. In Fig- 
ure|H| we show events in the energy interval 10 20 -^10 21 eV 
and in a ~ 0.1 ster solid angle, projected for three years 
of EUSO observations, assuming the apparent post-GZK 
AGASA flux (assumed to be isotropic). This is to be 
compared with the clustering-enhanced Z-burst emis- 
sions from the Virgo cluster, also shown. The latter con- 
stitute up to 3 % and 20 % of the projected primary 
events for m Vi — 1.5 eV and 0.6 eV, respectively, and 
dominate, albeit at a lower rate, beyond 10 20 5 eV. Thus, 
detecting Z-burst emissions with an observatory larger 
than EUSO, even in the face of less than optimal UHECz; 
fluxes and a possible post-GZK cosmic ray background, 
may not be entirely hopeless if the Cz^B clusters appre- 
ciably. 



VI. CONCLUSION 

Resonant annihilation of ultra-high energy cosmic neu- 
trinos (UHECz^) on the cosmic neutrino background 
(CzTi) into Z bosons — the Z-burst mechanism — is a 
unique, albeit indirect, process capable of revealing the 
CiTi in its present state. The annihilation can be de- 
tected as absorption dips in the incident UHECz/ flux 
and/or as emission features in the Z-decay products 
(nucleons and photons) at energies above the Greisen- 
Zatsepin-Kuzmin (GZK) cut-off. General Z-burst ab- 
sorption and emission phenomenology, including the pos- 
sibility that Z-decay products may constitute the post- 
GZK ultra-high energy cosmic ray (UHECR) events ob- 
served by AGASA, has been investigated in a number of 
recent publications. In the present work, we have consid- 
ered exclusively the implications of a non-uniform CzaB 
for the Z-burst emission rates and sky map. In partic- 
ular, we have focussed on the effects of augmentations 
to the CzaB number density in and around large galaxy 
clusters due to gravitational clustering, and the poten- 
tial of future cosmic ray experiments, especially the Ex- 
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FIG. 8: Number of Z-burst events, for three quasi-degenerate 
neutrinos, originating within 9 — 10° from M87 expected at 
EUSO in three years, assuming the cascade limited UHECz^ 
fluxes. Triangles denote the nucleon events, while the squares, 
pentagons and hexagons refer respectively to the photon 
events for the minimal, moderate and high URB scenarios. 
The stars indicate the number of events observed by AGASA 



in the same energy range (10 



-1(T 



eV) projected for three 



years of EUSO observations in the same solid angle. 



treme Universe Space Observatory (EUSO), to observe 
clustering-enhanced Z-burst emission rates originating 
from the nearby Virgo cluster. 

The GZK suppression of nucleons from extragalactic 
sources provides a clean environment for the study of 
Z-burst initiated events at energies above ~ 10 20 eV. In- 
deed, for UHECi' fluxes saturating the electromagnetic 
cascade limit, Z-burst nucleon and photon fluxes origi- 
nating from a uniform Cz/B already dominate at these 
energies. Gravitational Cz^B clustering at the Virgo clus- 
ter further enhances the fluxes by up to a factor of sev- 
eral to more than forty, depending on the neutrino mass, 
at the centre of the cluster (Figure 0J . The enhance- 
ment decreases with angular distance from the centre in 
a manner highly sensitive to the neutrino mass (Figure 
[5J . This angular dependence of the emission events can 
in principle serve as an independent probe of the neutrino 



12 



mass, the shape of the Virgo halo, as well as the Z -burst 
mechanism as a whole. 

However, the statistics of Z-burst observation is neces- 
sarily limited by the magnitude of the available UHECV 
flux. For three years of EUSO observations, we find that 
a detection of clustering-enhanced Z-burst rates from 
the Virgo cluster is probable, provided that the UHECV 
fluxes are close to current observational limits (Figure 
EJ). Sources capable of generating such large fluxes will 
most likely involve physics beyond the standard model. 
A quasi-degenerate neutrino mass spectrum is also desir- 
able, since gravitational clustering is considerably more 
efficient for "large" neutrino masses, m Vi > 0.1 eV (Fig- 
ure IS}. Under these favourable conditions, one would 
also expect the associated absorption dips in the inci- 
dent UHECi' flux to be resolved by EUSO and/or other 
forthcoming dedicated UHECV detectors. 

In the case of less than optimal UHECz^ fluxes, e.g., 
fluxes below the electromagnetic cascade limit, or a hi- 
erarchical neutrino mass spectrum, the observation of 
enhanced Z-burst rates in the direction of Virgo seems 



to require an experimental sensitivity exceeding that of 
EUSO by at least two orders of magnitude (Figure 0. 
The Westerbork radio observatory in the Netherlands of- 
fers a tremendous new reach in UHEC^ and UHECR 
detection by looking for signals of their interactions with 
the lunar regolith. If it establishes a UHEC^ flux at en- 
ergies above 10 21 eV, then the final uncertainty in the 
-Z-burst mechanism is removed. This allows definite pre- 
dictions to be made for associated Z-burst emissions from 
the Virgo cluster. Such measurements would finally es- 
tablish the existence (or not) of the cosmic neutrino back- 
ground. 



Acknowledgments 

We thank Markus Ahlers and Birgit Eberle for useful 
discussions. T.J.W. acknowledges the support of the U.S. 
Department of Energy grant DE-FG05-85ER40226 and 
the NASA- ATP grant 02-0000-0151. 



[1] S. Hannestad, Ne w J. Phys. 6, 108 (2004) 
arXiv:hep-ph/0404239 . 

[2] T. J. Wcilcr. Phys. Rev. Lett. 49, 234 (1982); Astrophys. 
J. 285, 495 (1984). 

[3] E. Roulet, Phys. Rev. D 47, 5247 (1993); S. Yoshida, 
H. Y. Dai, C. C. Jui and P. Sommers, Astrophys. J. 
479, 547 (1997) [arXiv:astro-p h/9608186 ; B. Eberle, 
A. Ringwald, L. Song and T. J. Weiler, Phys. Rev. D 70, 
023007 (2004) arXiv:hep-ph/0401203 ; G. Barenboim, 
O. Mena Requejo and C . Quigg, ibld~71, 083002 (2005) 
arXiv:hep-ph/0412122 . 

[4] D. Fargion* B. Meleand A. Salis, Astrophys. J. 517, 
725 (1999) [arXiv:astro-ph7 9710029 ; T. J. Weiler, As - 
tropart. Phys. 11, 303 (1999) arXiv:hep-ph/9710431 ; 
S. Yoshida, G. Sigl and S. J. Lee, Phys. Rev. Lett. 
81, 5505 (1998) arXiv:hep-ph/9808324 ; Z. Fodor, 
S. D. Katz and A. R ingwald, ibid. 88, 171101 

(2002) |arXiv:hep-ph /0105064 ; JHEP 0206, 046 (2002) 
|arXiv:nep-ph/ 0203198| . 

[5] M. Takeda et al, Phys . Rev. Lett. 81, 1163 (1998) 
|arXiv: astro-ph /98071931 ; Ast ropart . Phys. 19, 447 

(2003) arXiv:astro-ph/0209422 . 

[6] D. J. Bird et al. [HiRes Collaboration], Astrophys. J. 424, 
491 (1994); 441, 144 (1995). 

[7] R. U. Abbasi et al. [High Resolution Fly's Eye 
Collaboration!, Phys. Rev. Lett. 92, 151101 (2004) 
arXiv:astro-ph/0208243 : T. Abu-Zayyad et al. [HiRes 
Collaboration] , arXiv:astro-ph/0208301 

[8] N. N. Efimov et al., in Proceedings of the astrophysical as- 
pects of the most energetic cosmic rays, World Scientific, 
Singapore, 1991, p. 20. 

[9] M. A. Lawrence, R. J. O. Reid and A. A. Watson, 
J. Phys. G 17, 733 (1991); M. Ave, J. A. Hinton, 
R. A. Vazquez, A. A. Watson and E. Zas, Phys. Rev. 
Lett. 85, 2244 (2000) |arXiv:astro- ph/0007386 . 
[10] K. Greisen, Phys. Rev. Lett. 16, 748 (1966); G. T. Zat- 
sepin and V. A. Kuzmin, JETP Lett. 4, 78 (1966) [Pisma 



Zh. Eksp. Teor. Fiz. 4, 114 (1966)]. 

[11] O. E. Kalashev, V. A. Kuzmin, D. V. Semikoz 
and G. Sigl, Phys Rev. D 65, 103003 (2002) 
arXiv:hep-ph/0112351 ; D. S. Gorbunov, P. G. Tinyakov 
and" S. V. Troitsky, Astrop art. Phys. 18, 643 
(2003) [arXiv:astro-ph/0206385l; D. V. Semikoz and 
G. Sigl, JCAP 0404, 003 (2004) arXiv:hep-ph/0309328 ; 
G. Gelmini, G. Varieschi and T. J. Weiler, Phys. Rev. D 
70, 113005 (20 04) |arXiv:hep -ph/0404272 . 

[12] I. Kravchenko, arXiv:astro-ph/0306408 

[13] P. W. Gorham, GL~Hebert, K. M. Liewer, C. J. Naudet, 
D. Saltzberg and D. Williams, Phys. Rev. Lett. 93, 
041101 (2004) |arXiv:astro-p h/0310232]. 

[14] R. M. Baltrusaitis et al., Phys. Rev. D 31, 2192 
(1985); S. Yoshida et al. [AGASA Collaboration], in Pro- 
ceedings of the 27th International Cosmic Ray Confer- 
ence, Hamburg, Germany, 2001,edited by K.-H. Kam- 
pert, G. Heinzelmann, and C. Spiering (Copernicus- 
Gesellschaft e.V.,Katlenburg-Lindau, 2001), Vol. 3, pp. 
1142-1145; L. A. Anchordoqui, J. L. Feng, H. Goldberg 
and A. D. Shapere, P hys. Rev. D 66, 103002 (2002) 
[arXiv:hep -ph /0207139 . 

[15] N. G. Lehtinen, P. W. Gorham, A. R. Jacobson and 
R. A. Roussel-Dupre, Phys. Rev. D 69. 013008 (2004) 
[arXiv:astro-ph/0309656| . 

[16] M. Ackermann et al., Astropart. Phys. 22. 339 (2005). 

[17] X. Bertou, P. Billoir, O. Deligny, C. Lachaud and 
A. Letessier-Selvon, Ast ropart. Phys. 17, 183 (2002) 
[arXiv:astro-ph /0104452] . 

[18] P. Gorham et al, [ANITA collaboration], Proposal 
SMEX03-0004-0019, April 19, 2003. 

[19] S. Bottai and S. Giurgola [EUSO Collaboration], in Pro- 
ceedings of the 28th International Cosmic Ray Confer- 
ence, Tsukuba, Japan, 2003,edited by T. Kajita, Y. 
Asaoka, A. Kawachi, Y. Matsubara, and M. Sasaki (Uni- 
versal Academy Press, Tokyo, 2003), pp. 1113-1116; 
S. Bottai [EUSO Collaboration], to appear in Proceed- 



13 



ings Incontro Nazionale di Astrofisica delle Alte Energie, [41] 

Roma, May 1 5-16, 2003. 

D. Saltzberg, |arXiv:astro-ph/0501364| 
P. Gorham, private communication. 

S. Barwick et al. [ANITA collaboration], Phys. Rev., in 
preparation. 

J. Bacelar, Talk given at the International ARENA 

Workshop, May 17 19, 2005, Zeuthen, Germany, 

http:/ /www-zeuthen. desy.de/nuastro/arena/ [42] 

Orbiting Wide-angle Light-collectors, http://owl.gsfc 

nasa.gov/ [43] 

P. Sreekumar et al [EGRET Collaboration], Astrophys. 

J. 494, 523 (1998) arXiv:astro-ph/9709257 . 

K. Mannheim, R. J. Protheroe and J. P. Rachen, Phys. 

Rev. D 63, 023003 (2001) |arXiv:astro-"p h/9812398]. [44] 

See Z. Fodor, S. D. Katz and A. Ringwald, JHEP in Q. [45] 

G. L. Fogli, E. Lisi, A. Marrone, D. Montanino, A. 
Palazz o and A. M. Rotunno, e Conf C030626, THAT05 
(2003) arXiv:hep-ph/0310012 . 

L. N. da Costa, W. Freudling, G. Wegner, R. Giovanelli, 

M. P. Haynes and J. J. Salzer, Astrophys. J. Lett. 468 

L5, (1996) arXiv:astro-ph/9606144 ; A. Dekel et al, As- [46] 

trophys. J. 522 1, (1999) arXiv:astro-ph/9812197 . 

A. Ringwald and Y. Y. Y. Wong, JCAP 0412, 005 (2004) 

arXiv:hep-ph/0408241 . [47] 
M. Tegmark et al. [SDSS Collaboration], Phys. Rev. D 
69, 103501 (2004) arXiv:astro-ph/0310723 . [48] 
M. Tegmark et al. [SDSS Collaboration], Astrophys. J. 
606, 702 (2004) |arXiv:a stro-ph /03107251 . [49] 
C. L. Bennett et al., Ast rophys. J. Suppl. 148, 1 (2003) 

arXiv:astro-ph/0302207 . 

(J. Weinheimer, JNucl. Phys. Proc. Suppl. 118, 279 [50] 
(2003); V. M. Lobashev et al, ibid. 91, 280 (2001). 

H. V. Klapdor-Kleingrothaus et al, Eur. Phys. J. A 12, 

147 (2001) arXiv:hep-ph/0103062 ; C. E. Aalseth et al. [51] 
[IGEX Collaboration], Phys. Rev. D 65, 092007 (2002) 
arXiv:hep-ex/0202026 . 

O. Elgaro y and O. Lah av, New J. Phys. 7, 61 (2005) [52] 
arXiv:hep-ph/0412075 . 

S. Singh and C P. Ma, P hys. Rev. D 67, 023506 (2003) [53] 
arXiv:astro-ph/0208419 . [54] 
P. Fouque et al, Astron . Astrophys. 375, 770 (2001) [55] 
|arXiv:astro-ph/010626"T1. [56] 
http:/ /www. euso-mission.org/ 
http://www.auger.org/ 



C. Lunardini and A. Yu. Smirnov, Phys. Rev. D 64, 
073006 (2001) |arXiv:hep-ph/0012056| ; A. D. Dolgov, 
S. H. Hansen, S. Pastor, S. T. Petcov, G. G. Raf- 
felt and D. V Semikoz ,Nucl. Phys. B 632, 363 (2002) 
|arXiv:hep-p h/02012 87|; Y. Y. Y. Wong Phys. Rev. D 
66, 025015 (2002) arXiv:hep-ph/0203180 ; K. N. Abaza- 
jian, J. F. Beacom and N F. Bell, ibid. 66, 013008 (2002) 
[arXiv:astro -ph /0203442 . 

A. D. Dolgov and F. Takahashi, Nucl. Phys. B 688, 189 
(2004) arXiv:hep-ph/0402066 . 

Z. Fod or and S. D. KatzTPh ys. Rev. D 63, 023002 
(2001) arXiv:hep-ph/0007158 ; Z. Fodor, S. D. Katz, 
A. Ringwald and H. Tu, Phys. Lett. B 561. 191 (2003) 
arXiv:hep-ph/0303080 . 
http://www.desy.de/~uhecr 
S. Lee, Phys. R~ D 
arXiv:astro-ph/9604098 ; R. 

A. Johnson, Astropart. 
arXiv:astro-ph/9506119 ; R. J 



58, 043004 (1998) 
J. Protheroe and 
Phys. 4, 253 (1996) 
Protheroe and P. L. Bier- 



mann, ibid. 6, 45 (1 996) [Erratum, ibid. 7, 181 (1997)] 
arXiv:astro-ph/9605119 . 

j. F. Navarro, C. S. Fre nk and S. D. M. White, A stro- 
phys. J. 462, 563 (1996) |arXiv:astro-ph /9508025 ; 490, 
493 (1997). 

J. S. Bullo ck et al, Mon. Not. Roy . Astron. Soc. 321, 
559 (2001) |arXiv:astro-p h/9908 1 59] . 
V. S. Berezinsky and A. Yu. Smirnov, Astrophysics and 
Space Science 32, 461 (1975). 

U. Kes het, E. Waxman and A. Loeb, JCAP 0404, 006 

(2004) [arXiv:astro-'p h/030644 2]; A. W. Strong, I V. 
Moskalenko and O. Reimer, arXiv:astro-ph/0306345 

M. Ahlers, L. A. Anchordoqui, H. Goldberg, F. Halzen, 
A. Ringwald an d T. J. Weiler, Phys. Rev. D 72, 023001 

(2005) arXiv:astro-ph/0503229 . 

V. Berezinsky, in Texas in Tuscany, Proceedings of the 
21st Texas Symposium on relativistic astrophysics, World 
Scientific, River Edge NJ, 2003 arXiv:hep-ph/0303091 . 
V. S. Berezinsky and A . Vilenkin, P hys. Rev. D 62, 
083512 (2000) arXiv:hep-ph/9908257 . 
See O. E. Kalashev et al. in ITU . 
R. Gandhi et al, Phys. Rev. D 58, 093009 (1998). 
See B. Eberle et al in Q. 

H. Falcke et al. [LOPES Collaboration], Nature 435, 313 
(2005); see also |http: / /www.astro . ru.nl/lopes / 



